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Abstract 

In this review article we discuss recent constructions of global F-thcory GUT models and 
explain how to make use of toric geometry to do calculations within this framework. After 
introducing the basic properties of global F-theory GUTs we give a self-contained review 
of toric geometry and introduce all the tools that are necessary to construct and analyze 
global F-theory models. We will explain how to systematically obtain a large class of 
compact Calabi-Yau fourfolds which can support F-theory GUTs by using the software 
package PALP. 
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1 Introduction 

Even though it has been around for quite a while PQ, F-theory has recently received a lot 
of new attention as a setup where Grand Unified Theories (GUTs) can be conceived from 
string theory. Starting with [21 [3l d] the phenomenology of F-theory GUTs has become an 
active field of research. The basic idea is that the GUT theory is localized on a (p, q) seven- 
brane S inside a three-dimensional base manifold B of an F-theory compactification on an 
elliptically fibered Calabi-Yau fourfold. The location of the GUT brane and the gauge group 
are determined by the degeneration of the elliptic fibration. Chiral matter localizes on curves 
inside the GUT brane S where gauge enhancement occurs, Yukawa couplings sit at points. 
For many phenomenological applications it is sufficient to consider the field theory living 
on the GUT brane without specifying the details of the global F-theory compactification. 
However, fluxes, monodromies or consistency constraints such as tadpole cancellation cannot 
be addressed in a purely local setup. These issues have recently received a lot of attention 
in the literature El El O El [TOl [IH H21 [13 [TH [151 H3 HZl [13 H21 [20l EJJ [22l [23l EH E5]. 
Therefore it is interesting to see whether it is possible to embed the local F-theory GUT 
into a compactification on a Calabi-Yau fourfold. Most known examples of compact Calabi- 
Yau manifolds are hypersurfaces or complete intersections in a toric ambient space. It is 
thus natural to look for Calabi-Yau fourfolds within this class of examples. A prescription 
for constructing elliptically fibered Calabi-Yau fourfolds as complete intersections in a six- 
dimensional toric ambient space has been given in [SI IS]- Before that complete intersection 
Calabi-Yau fourfolds in F-theory had already been used in the context of F-theory uplifts of 
type IIB string theory [26} |2"71 128] . A similar construction has also been discussed in |l(Jj . 
It has been shown in examples that it is indeed possible to construct viable F-theory GUTs 
within this framework. 

The construction of [6] is very well-suited for a systematic search of a large class of models. 
This is interesting for several reasons: one goal is to find particularly nice examples of F- 
theory compactifications. Even though the known examples have been able to incorporate 
F-theory models, one usually gets much more than just that. In minimal F-theory GUTs one 
typically needs only very few Yukawa points and a small number of moduli on the matter 
curves. This is not satisfied in most known global models. A related question deals with the 
genericity of F-theory GUTs. The geometric configurations used for constructing such models 
are usually quite special and one may wonder how often they can be realized in elliptically 
fibered fourfolds. From the point of view of model building it is useful to have some easy- 
to-check geometric conditions which makes it possible to select suitable models from a large 
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class of geometries. This will be discussed in more details in the text. From a mathematical 
point of view it might be interesting to obtain a partial classification of Calabi-Yau fourfolds. 

This review article discusses selected topics in toric geometry and F-theory GUTs. The 
paper is organized as follows: in section [2] we recall the construction of global F-theory models 
and discuss the basic requirements we would like to impose. In section [3] we review several 
notions in toric geometry which are required in order to perform the F-theory calculations. 
The geometries one has to deal with are usually quite complicated, and very often one has 
to rely on computer support in order to be able to do explicit calculations. Therefore we 
discuss how such calculations can be implemented using existing software such as PALP [29]. 
We will mainly focus on the application of toric geometry in the context of F-theory model 
building. For a more complete picture on this vast subject of F-theory phenomenology we 
refer to other review articles such as [301 Ell E2]- For more extensive discussions of toric 
geometry we recommend [331 EH E5] • 

2 Global F-theory models 
2.1 Setup 

In this section we introduce the basic concepts and notions used in global F-theory models. 
In the remainder of this review we will explain the techniques which are necessary to do 
calculations within this framework. For more details on how the quantities introduced below 
come about we refer to the original papers or the recent review article [32] . 

In [B] it has been proposed to construct Calabi-Yau fourfolds, which are suitable for F- 
Theory model building, as complete intersections of two hypersurfaces in a six-dimensional 
toric ambient space. The hypersurface equations have the following structure: 

PB(yi, w) = Pw(x, y, z, y h w) = (1) 

The first equation only depends on the homogeneous coordinates (?/,, w) of the three-dimensional 
base B of the elliptically fibered Calabi-Yau fourfold X4. Here we have singled out one coor- 
dinate w, indicating that the divisor given by w = defines a seven-brane S which supports 
a GUT theory of the type introduced in [2 [31 H] . The second equation in ([1]) defines a Weier- 
strass model, where (x, y, z) are those coordinates of the six-dimensional ambient space that 
describe the torus fiber. For this type of elliptic fibrations Py/ is of Tate form which is defined 
as follows: 

P\y = x 3 — y 2 + xyzai + x 2 z 2 a2 + yz 3 a^ + xz^a^ + z 6 ciq (2) 

The a n (yi,w) are sections of K^ n , where Kb is the canonical bundle of the base manifold. 
Furthermore x and y can be viewed as sections of Kg 2 and Kg 3 , respectively. 

The information about the F-theory model is encoded in the Tate equation ([2]). In order 
to have a non-trivial gauge group on the GUT brane the elliptic fibration must degenerate 
above S. The gauge group is determined by the structure of the singularity. The elliptic 
fibration becomes singular at the zero locus of the discriminant A. Defining the polynomials 
P2 = of + 4a2, Pi = 0103 + 2a4 and {3$ = a 2 + Aa% the discriminant is given by the following 
expression: 

A = - 01) - 8/3f - 27$ + 9/3 2 /3 4 /3 6 (3) 
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According to Kodaira's classification [36] and Tate's algorithm |37j . the gauge group can be 
inferred from the factorization properties of the a n (yi,w) with respect to w. Considering for 
instance SU(5)- and 5O(10)-models, the factorization looks like this: 

SU (5) : ai = 65^° C12 = b^w 1 03 = b^w 2 04 = 62W 3 a§ = bow 5 
50(10) : a\ = b^w 1 a>2 = b^w 1 03 = b^w 2 04 = 62W 3 a,Q = bow 5 (4) 

The bis are sections of some appropriate line bundle over B that have at least one term 
independent of w. 

In F-theory GUT models chiral matter localizes on curves on 5, where a rank 1 enhance- 
ment of the gauge group appears. In SU(5) models the matter curves are at the following 
loci inside 5: 

6364 — 626365 + bob 2 = 5 matter 65 = 10 matter , (5) 

The matter curves for the 50(10) models are at: 

63 = 10 matter 64 = 16 matter . (6) 

Yukawa coupling arise at points inside 5 where the GUT singularity has an rank 2 enhance- 
ment. In SU(5) models the Yukawa points sit at: 

64 = 01165 = 10105 Yukawas Eq enhancement 

65 = n 63 = 10 5 5 Yukawas 50(12) enhancement . (7) 

In the 50 (lO)-case we have the following Yukawa couplings: 

63 = 0(164 = 161610 Yukawas £7 enhancement 

6| - 46 6 4 = n 63 = 16 10 10 Yukawas 50(14) enhancement . (8) 

Given a complete intersection Calabi-Yau fourfold of the form ([1]) the expressions for matter 
curves and Yukawa points are globally defined and can be calculated explicitly. Having a 
global F-theory compactification we can furthermore calculate the Hodge numbers and the 
Euler number X4 of the Calabi-Yau fourfold X4. The latter enters the D3-tadpole cancellation 
condition, 

g = iV D 3 + ^G 4 AG 4 , (9) 
where G4 denotes the fourform flux on X4 and Nps is the number of D3-branes. 

2.2 Geometric Data in F-theory models 

So far we have summarized the basic structure of a global F-theory GUT. In the present 
section we will discuss which properties of the GUT model are encoded in the geometries 
of the base manifold B and the Calabi-Yau fourfold X4. We will not go deeply into the 
phenomenology of F-theory GUTs but rather focus on the basic geometric properties which 
should be satisfied in order to obtain a viable GUT model. 
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Base Manifold 



Since the GUT brane S is a divisor in a three-dimensional base manifold B, a large amount 
of information about the model can be extracted from the geometry of B. The base B is a 
non-negatively curved manifold of complex dimension three. In our setup it will be given by 
a hypersurface in a toric ambient space. Note that Fano threefolds are not a good choice for 
B due to the lack of a decoupling limit [38] . In section [3] we discuss a systematic construction 
of such base manifolds using toric geometry. In order to have a well defined model we have 
to make sure that B is non-singular. In contrast to Calabi-Yau threefolds the base manifolds 
for F-theory GUTs may inherit the singularities of the ambient space. Therefore checks for 
the regularity of B have to be implemented. 

Having found a suitable base manifold the next step is to identify divisors inside B that 
can support GUT models. The most promising candidates for F-theory model building are 
del Pezzo surfaces. These are Fano twofolds (see for instance [39]). Note however, that del 
Pezzos are not the only possibility for the construction of GUT models in F-theory. See [ID] 
for a recent discussion. There are several motivations to focus on del Pezzo divisors. In local 
F-theory GUTs the del Pezzo property ensures the existence of a decoupling limit OH]. For 
577(5) GUT models, the fact that del Pezzos have h ' 1 = h 2 ' = implies some powerful 
vanishing theorems which forbid exotic matter after breaking 577 (5) to the Standard Model 
gauge group [I]. 

We can identify candidates for del Pezzo divisors inside B by their topological data. 
Suppose the base manifold B is embedded in a toric ambient space which has toric divisors 
Dj. The Di give a homology basis of the ambient space. In this setup the hypersurface is 
specified by a divisor, which we will by abuse of notation also call B, that is given in terms of 
a linear combination of the Di. The total Chern class of a particular divisor 5 in the ambient 
space is, after restriction to B (for more details see section [33]) : 

c(S ) - n,(i+A) , m 

In order to apply this formula we have to know the intersection ring of B. As we will discuss 
in section [3731 this can be obtained from the intersection ring of the ambient space. 

A necessary condition for a divisor 5 to be dP n is that it must have the following topo- 
logical data: 

[ Cl (S) 2 = 9-n [ C2 (S) = n + 3 X h = / Td(5) = 1, (11) 

Js Js Js 

where \h = 1)*^°'* is the holomorphic Euler characteristic and Td denotes the Todd 

class. In the equations above the integration over the four-cycle (representing the divisor) 5 
is equivalent to computing the intersection with S. Since del Pezzos are Fano twofolds, we 
have a further necessary condition: the intersections of c\(S) with curves on 5 have to be 
positive. In the toric setup we can only check this for curves which are induced from the 
divisors on the ambient space. In that case the condition is: 

Di ■ S ■ ci(5) > Di ^5 VA • S i- . (12) 

In order to make these calculations we need to know the homology basis of toric divisors and 
their intersection numbers. 
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In local F-theory GUTs the del Pezzo property is sufficient to ensure the existence of a 
decoupling limit. For global models further checks are in order. Gravity decouples from the 
gauge degrees of freedom if the mass ratio Mqut/M p i becomes parametrically small. The 
Planck mass M p \ and the mass scale Mqut of the GUT theory are related to the geometry 
of B and S in the following way: 

M 8 i M 4 
A$~-fVolCB) M gut ~ Vol(S)-* 1/<7 YM 2-Vol(5), (13) 



9% 9s 



Therefore one has: 



M GUT 2 Vol(S) 3 / 4 
"M^^Vol^ ^ 

There are two ways to achieve a small value for Mgut/M p i, now commonly referred to as the 
physical and the mathematical decoupling limit. In the physical decoupling limit the volume 
of the GUT brane S is kept finite while Vol(-B) — > oo. The mathematical decoupling limit 
takes Vol(S') — > for finite volume of B. The two limits may not be equivalent in the sense 
that they may be implemented by tuning different Kahler parameters. The volumes of B and 
S can be determined in terms of the Kahler form J of the ambient toric variety restricted to 
B: 

Vol(-B) = J 3 Vol(S) = S ■ J 2 (15) 

In order for the volumes to be positive we must find a basis of the Kahler cone Ki where, by 
definition, J can be written as J = Yli r iKi with n > 0. The existence of mathematical and 
physical decoupling limits can be deduced from the moduli dependence of these volumes. 

Having found a suitable base manifold we can also study matter curves and Yukawa 
couplings. The curve classes M of the matter curves can be expressed in terms of the toric 
divisors of the ambient space. The genus of the matter curves can be computed using the 
first Chern class of the matter curve and the triple intersection numbers: 

C(M) = (1 + B)(1 + 5)(1 + M) (16) 

Here we have assumed that M is irreducible. After expanding this expression to get ci(M), 
the Euler number is obtained by the following intersection product: 

X (M) = 2- 2g(M) = ci(M) ■ M ■ S (17) 

The genus of a matter curve gives us information about the number of moduli the curve has. 
Since these moduli have to be stabilized, matter curves of low genus are desirable from a 
phenomenological point of view. In the generic situation the equations specifying the Yukawa 
points can be expressed as classes Y\, Y% in terms of the toric divisors. The number of Yukawa 
points is then given by the following intersection product: 

"-Yukawa = S • Y\ ■ Y 2 (18) 

In order to account for the Standard Model Yukawa couplings only a small number of Yukawa 
points is needed. 
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Fourfold 

Given a base manifold B one can construct a Calabi-Yau fourfold X4 which is an elliptic 
fibration over B. As described in the next section this can be done systematically using toric 
geometry. However, not all of the desirable features of global F-theory models are automatic 
in this construction. The main requirement on X4 is that it is a complete intersection of 
two hyper surf aces. Furthermore these hypersurfaces must have a specific structure (P). In 
order to be able to use powerful mathematical tools we furthermore have to require that there 
exists a nef-partition (cf. section I3.2[) which is compatible with the elliptic fibration. When 
this elementary requirement is satisfied we can engineer a GUT model. This is done in two 
steps: first we have to identify the GUT divisor S, given by the equation w = in B within 
the Calabi-Yau fourfold. The second step is to impose the GUT group. This amounts to 
explicitly imposing the factorization conditions such as (HD on the Tate model. This means 
that we have to remove all those monomials in (J2|) which do not satisfy the factorization 
constraints. This amounts to fixing a number of complex structure moduli of X4. 

Recently there has been active discussion in the literature how to globally define fluxes in 
F-theory models [SI HH [151 ESI [201 EH EH]. In F-theory model building fluxes enter at 
several crucial points. Gauge flux along the matter curve is needed in order to generate chiral 
matter. Breaking of the GUT group to the Standard Model gauge group can be achieved by 
turning on ?7(l)-flux. Furthermore, in SU(5) F-theory GUTs we need global U(l)s in order 
to forbid dimension 4 proton decay operators. In SO (10) F-theory GUTs they are needed 
in order to obtain chiral matter |41| I13j . A general global description of the fourform flux 
G4 is still missing. In [12] an auxiliary construction involving spectral covers that factorize 
was used to describe fluxes locally in the vicinity of the GUT brane. It has been shown 
in |14t \15\ [16] that under certain circumstances the information captured by the spectral 
cover can be encoded in the Tate model, and is therefore global. However, this need not be 
the case |llj . In [2] it has been shown that a spectral cover which factorizes is generically 
globally defined for "t7(l)-restricted Tate models". This is achieved by imposing a global 
U(l)x symmetry in the elliptic fibration. In terms of the Tate model this is achieved by 
setting ci6 = 0. 

3 Ingredients and Techniques from Toric Geometry 

In the previous section we have introduced quantities which encode important information 
about F-theory GUT models in the geometry of the base manifold and the Calabi-Yau four- 
fold. In this section we will provide the tools to calculate them. The input data needed for 
these calculations can be obtained by using toric geometry. After giving the basic definitions 
we will discuss how to describe hypersurfaces and complete intersections of hypersurfaces in 
toric ambient spaces. Then we explain how to obtain the intersection ring and the Kahler 
cone, or dually, the Mori cone. Finally we will discuss how to use the computer program 
PALP [29] for calculations in toric geometry. This discussion of toric geometry has been 
compiled with a view towards the applications in F-theory model building. It is by no means 
an exhaustive description of this vast subject which brings together algebraic geometry and 
combinatorics. 
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3.1 Toric Varieties 



We start by defining a toric variety X of dimension n as the following quotient: 

X = (C r - Z)/((C*) r - n x G), (19) 

where G is a finite abelian group, (C*) r_n describes the action of an algebraic (r — n)-torus 
and Z C C r is an exceptional set which tells us which combinations of coordinates are not 
allowed to vanish simultaneously. The simplest example is CP 2 , where the C*-action is given 
by (zi,Z2, zs) ~ (Xz±,Xz2,Xzs), A G C*, the exceptional set is Z = {z\ = Z2 = ^3 = 0} 
and G is trivial. Thus, as is well-known: CP 2 = (C 3 — {z± = Z2 = z% = 0})/( (z\, £2,23) ~ 
(Azi, Az 2 ,Az 3 )). 

The crucial fact about toric geometry is that the geometric data of the toric variety can 
be described in terms of combinatorics of cones and polytopes in dual pairs of integer lattices. 
The information about the toric variety is encoded in a fan E, which is a collection of strongly 
convex rational polyhedral cones where all the faces and intersections of pairs of cones also 
belong to the fan. 'Strongly convex' means that all the cones of the fan have an apex at 
the origin, 'rational' means that the rays that span the cone go through points of the lattice. 
We denote by £( n ) the set of all n-dimensional cones. In order to define the fan we use the 
fact that a toric variety X contains an n-torus T = (C*) n as a dense open subset whose 
action extends to X. Parametrizing T by coordinates (t\, . . . ,t n ), one defines the character 
group M = {x : T — > C*} and the one-parameter subgroups N = {A : C* — > T}. M and 
N can be identified with integer lattices that are isomorphic to TU l . Given a point m G M, 
the character is given by x m {t) = t™ 1 ■ ■ - t™ n = t m . This is a holomorphic function on T, 
and descends to a rational function on the toric variety X. For every u G N, A is defined as 
A"(t) = (t Ui , . . . , T Un ) for r G C*. The fan S and its cones a are defined on the real extension 
-ZVr of N. The lattices M, N are dual due to the composition (x o A)(r) = x(A(r)) = t^ x ' X \ 
where (x m , A") = m ■ u is the scalar product. 

The M-lattice encodes the data about regular monomials in X, the iV-lattice captures 
the information about the divisors. The divisors defined by x m = can be decomposed in 
terms of irreducible divisors Dj\ div(x' m ) = X^=i a 3^y These divisors are principal divisors, 
i.e. divisors of meromorphic function where Da correspond to poles or zeros and the aj are 
orders of the pole/zero. The coefficients aj(m) G Z are unique, and there exists a map 
m — > a,j(m) = (m,Vj) with Vj G N. Thus there is a vector Vj for every irreducible divisor 
Dj. The Vj are the primitive generators of the one-dimensional cones pj (i.e. rays) in the fan 
S. The convex hull of the Vj defines a polytope A* = conv{vj}. Locally, we can write the 
divisors as Dj = {zj = 0}, where Zj is regarded as a local section of a line bundle. Dj are 
called toric divisors. There are linear relations among the Vj G A* which translate into linear 
relations among the toric divisors. 

In order to make contact with the definition (|19p of X, we view the {zj} as global ho- 
mogeneous coordinates (z\ : ... : z r ). If all Zj are non-zero the coordinates (A 9l zi : ... : 
X qr z r ) ~ (z\ : . . . : z r ) with A G C* describe the same point on the torus T, if Yl ( lj v j = 
for Vj G N as above. Since the Vj live in an n-dimensional lattice they satisfy r — n linear 
relations. If the Vj do not span the iV-lattice there is a finite abelian group G such that 
G ~ N/(sp&n{vi, . . . ,v r }). Identifications coming from the action of G have to be added 
to the identifications between the homogeneous coordinates coming from the torus action. 
Having introduced the fan S, we are also able to specify the exceptional set Z that tells 
us where the homogeneous coordinates are not allowed to vanish simultaneously: a subset 
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of coordinates Zj is allowed to vanish simultaneously if and only if there is a cone a G E 
containing all the corresponding rays pj. To be more precise, the exceptional set is the union 
of sets Zj with minimal index sets / of rays for which there is no cone that contains them: 
Z = UiZj. This is equivalent to the statement that the corresponding divisors Dj intersect 
in X. Putting the pieces together we arrive at the definition (|19j) . 

There are two important properties of the fan E which translate into crucial properties 
of the toric variety X. Firstly, X is compact if and only if the fan is complete, i.e. if the 
support of the fan covers the iV-lattice: |£| = |J S a = Nr. Secondly, X is non-singular if and 
only if all cones are simplicial and basic, which means that all cones a G E are generated by 
a subset of a lattice basis of N. Singularities can be removed by blow-ups, where singular 
points are replaced by P n_1 s. All the singularities of a toric variety can be resolved by a series 
of blowups. These correspond to subdivisions of the fan. In order to completely resolve all 
singularities one must find a maximal triangulation of the fan. In many cases it is sufficient 
to find a maximal triangulation of the polytope A*. 

Finally, let us emphasize the significance of the homogeneous weights q{. In general there 
will be a full (r — n) x r matrix Qij, called weight matrix, whose (r — n) lines encode the 
C*-actions. Since each of the Zj corresponds to an irreducible divisor in X, the columns of 
the weight matrix define a homology basis of the divisors Dj . In physics language the weights 
qi are the £7(l)-charges in the gauged linear sigma model that defines the toric variety X. 
Note that the weights contain all the information to recover the M- and iV-lattice. With the 
weight matrix as input this can be done using PALP. 

3.2 Hypersurfaces and complete intersections 

Having defined a toric variety we go on to discuss hypersurfaces and complete intersections 
of hypersurfaces in toric varieties. The hypersurface equations are sections of non-trivial line 
bundles. The information of these bundles can be recovered from their transition functions. 
In this context we introduce the notions of Cartier divisors and Weil divisors. A Cartier 
divisor is given, by definition, by rational equations f a = and regular transition functions 
fa/f/3 on the overlap of two coordinate patches U a ,Up. Cartier divisor classes determine the 
Picard group Pic(A) of holomorphic line bundles. Weil divisors are finite formal sums of 
irreducible varieties of codimension 1. On a toric variety the Chow group A n _\{X) modulo 
linear equivalence is generated by the T-invariant irreducible divisors Dj modulo the principal 
divisors div(x m ), m G M. A Weil divisor of the form D = ^CLjDj is Cartier if there exists 
an m a G M for each maximal cone a G E( n ) such that (m a ,Vj) = —dj for all rays pj G a. If 
X is smooth then all Weil divisors are Cartier. If X is compact and D is Cartier then O(D) 
is generated by global sections if and only if (m a , Vj) > —a,- for a G £( n ) and pj <£_ a. If this 
is the case for v G a, iPd{v) — ( m ai v ) is a strongly convex support function. With that we 
can define a polytope Ajj = {m G Mr : (m a ,Vj) > — dj}. This is a convex lattice polytope 
in Mk whose lattice points provide global sections of the line bundle 0{D) corresponding to 
the divisor D. D is generated by global sections if and only if Ad is the convex hull of {m CT }. 
Furthermore, D is ample if and only if Ad is n-dimensional with vertices m CT for a G £(") and 
with m a ^ m T for a 7^ r G E^ n ). Finally D is called base point free if and only if m a G Ad for 
all a G £( n ). Base point freedom is a sufficient condition for a hypersurface defined by D to 
be regular: Bertini's theorem states that the singular points of D are the base locus and the 
singular points inherited from the ambient space. The absence of base points implies that D 
can be deformed transversally in every point and therefore generically avoids the singularities 
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of the ambient space. Thus, a base point free D is regular. We emphasize however that base 
point freedom is not a necessary condition for the regularity of D. 

Equations for hypersurfaces or complete intersections are sections of line bundles 0(D) 
given by the following Laurent polynomial: 

/= £ c mX m = £ c m n^ J> (20) 

rneA D nM m&A D nM j 

In an affine patch U a the local section f a = //x m<T is a regular function. Given a polytope 
Ad G M, we can define the polar polytope A° D by A° D = {y G iNfo : (x, y) > -1 Vi G Ad}. 
It can be shown [33] that the Calabi-Yau condition for hypersurfaces requires that Ad Q Mr 
is polar to A* = A° D C Ar, where A* is the convex hull of the Vj G N as defined in section 
13.11 A lattice polytope whose polar polytope is again a lattice polytope is called reflexive. 
For reflexive polytopes (A, A°) there exists a combinatorial formula for the Hodge numbers 
03]: 

hi s i(X A ) = h dimA _ 2A (X A o) = Z(A°) - 1 - dimA - £ l*(9°) + ]T l*(9°)l*(9) 

codim(6»°)=l codim(6»°)=2 

(21) 

where 9 and 9° is a dual pair of faces of A and A . Furthermore, 1(9) is the number of lattice 
points of a face 9 and I* (9) is the number of its interior lattice points. 

In our discussion of F-theory model building we also encounter complete intersection 
Calabi-Yaus. The concept of polar pairs of reflexive lattice polytopes can be generalized as 
follows: 

A = A! + ... + A r A° = (Vi,...,V r ) CO nv 

(V n) A m ) > -5 nm (22) 
V° = (Ai, . . . , A r ) conv V = Vi + . . . + V r 

Here r is the codimension of the Calabi-Yau and the defining equations fi = are sections 
of O(Ai). The decomposition of the M-lattice polytope A C Mr into a Minkowski su 
A = Ai + . . . + A r is dual to a nef (numerically effective) partition of the vertices of a 
reflexive polytope V C Ar such that the convex hulls (V;) conv of the respective vertices and 
G A only intersect at the origin. The nef-property means that the restriction of the line 
bundles associated to the divisors specified by the N-lattice points to any algebraic curve of 
the variety are non-negative. There exists a combinatorial formula for the Hodge numbers 
[44] which has been implemented in PALP. 

In many string theory applications, and in particular also in F-theory, the fibration struc- 
ture of a Calabi-Yau manifold is of great interest. For Calabi-Yaus which can be described 
in terms of toric geometry the fibration structure can be deduced from the geometry of the 
lattice polytopes. If we are looking for toric fibrations where the fibers are Calabi-Yaus of 
lower dimensions, we have to search for reflexive sub-polytopes of A° which have appropriate 
dimension. Given a base b and a fiber /, the fibrations descend from toric morphisms of the 
ambient spaces corresponding to a map (f) : £ — >■ of fans in and Nj,, where eft : N — >■ Nf, is 
a lattice homomorphism such that for each cone a G 5] there is a cone G Sf, that contains 
the image of a. The lattice Nf for the fiber is the kernel of 4> in N. The fiber polytope is then 

lr The Minkowski sum A + B of two sets A, B is defined as follows: A + B = {a + b\a G A, b G B}. 
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defined as follows: A^ = A° n Nf. In order to guarantee the existence of a projection one 
must find a triangulation of A j and extend it to a triangulation of A° . For each choice of tri- 
angulation the homogeneous coordinates corresponding to the rays in A^ can be interpreted 
as coordinates of the fiber. 

3.3 Intersection ring and Mori cone 

Two further pieces of data that are necessary in many string theory calculations are the 
intersection numbers of the toric divisors and the Mori cone, which is the dual of the Kahler 
cone. Inside the Kahler cone the volumes such as (|15p are positive. Thus, in the context 
of F-theory model building the Kahler cone is needed in order to make statements about a 
decoupling limit. 

Let us start with discussing the intersection ring. For a compact toric variety the 
intersection ring is of the form Z[Di, . . . , D r ]/{Iu n + I non -n n ). The two ideals to be divided 
out take into account linear and non-linear relations between the divisors. The linear relations 
have the form ^2j{m,Vj)Dj, where m £ M form a set of basis vectors in the M-lattice. The 
non-linear relations are denoted by R = URj where the Rj are of the form Rj = Dj 1 -. . ■■Dj k = 
0. They come from the exceptional set Z = L)Zj defined in section 13.11 which determines 
which homogeneous coordinates are not allowed to vanish at the same time. As mentioned 
before, this is the case when a collection of rays Pj t ,. ■ ■ , Pj k € N is not contained in a single 
cone. The non- linear relations R generate the ideal I non -iin which is called Stanley- Reisner 
ideal. Thus, the intersection ring A*(E) of a non-singular toric variety has the following form: 

A4Xx) = Z[D 1 ,...,D r ]/(R,J2(m,v j )D j ) (23) 

3 

The definition of the intersection ring holds for non-singular toric varieties but may be gen- 
eralized to the case where Xy, is simplicial projective. This means that the toric variety may 
be singular but still all the cones of the fan E are simplicial. Such a situation may occur for 
example if we choose a non-maximal triangulation of the polytope A*. In this case the inter- 
section numbers take values in Q. To compute the Stanley-Reisner ideal in the non-singular 
case one must find a maximal triangulation of the fan E or the polytope A*. In order to 
get intersection numbers we still have to fix a normalization: for a maximal simplicial cone 
a € E( n ) spanned by vj 1 , . . . , Vj n we fix the intersection numbers of the corresponding divisors 
to be Dj t ■ . . . Dj n = l/Vol(<r), where VoI(ct) is the lattice volume of a (i.e. the geometric 
volume divided by the volume 1/n! of a basic simplex). If X is non-singular the volume is 1. 
Using the intersection ring one can compute the total Chern class of the tangent bundle Tx 
of X which is given by the following formula: c(Tx) = Hj=i(l + Dj). 

So far, we have only discussed the intersection ring of the toric variety X. However in 
many applications we rather need the intersection numbers for divisors on a hypersurface 
given by a divisor D in X. Here we can make use of the restriction formula that relates the 
intersection form on the hypersurface divisor to the intersection form on X: 

D jn _, |d = Dj 1 ■ . . . ■ D jn _, ■ D\ x (24) 

This allows us to compute the intersection ring of D from the intersection ring of X. In 
([23]) restriction to D amounts to computing the ideal quotients of Ii in and I non -iin with the 
ideal generated by D. By adjunction the Chern class for the hypersurface specified by D is 
c{D)=YF j=1 {l + D j )/{l + D). 
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In order to be able to calculate all the quantities denned in section 12.21 we miss one more 
ingredient: the Mori cone. By definition, the Mori cone is the dual of the Kahler cone. We 
need the information about the Kahler cone in order to be able to compute the volumes 
of divisors. By definition the volumes will be positive inside the Kahler cone. The Mori 
cone is generated by l^ 1 where k = r — n if the fan E is simplicial. Otherwise 
the number of Mori generators can be larger. The Mori cone L is then defined as follows: 
+ . . . + M>o^ fc ). For the calculation of the Mori cone we also require a maximal 
triangulation of A*. Given such a triangulation the Mori generators can be determined as 
follows |45j : take every pair of n-dimensional simplices (Sk, Si) which have a common n — 1- 
dimensional simplex s^i = S^HSi. Then find the unique linear relation ^ l i ' Vi = with vi € 
SfrUSi where the l^' 1 are minimal integers and the coefficients of the points in (SkUSi)\(Sf.C\Si) 
are non- negative. The Mori generators are then the minimal integers by which every l k ' 1 
can be expressed as positive integer linear combinations. There is an equivalent algorithm to 
determine the Mori generators due to Oda and Park [36] which has been implemented in an 
unreleased version of PALP [37]. Note that the relations Yll=i = define the ideal Iu n 

in (|23p . Assembling the Mori vectors into a k x r-matrix, the columns of the matrix encode 
inequalities for the values of the Kahler parameters. Solving these inequalities yields a basis 
Ki of the Kahler cone such that the Kahler form of X can be written as J = riKi with 
r, > 0. Note that this prescription computes the Kahler cone of the toric variety X. It is 
often assumed that this is a good approximation for the Kahler cone of a hypersurface in X. 

3.4 Toric Calculations using PALP and other Software 

In string theory and F-theory we deal with compactifications on Calabi-Yau threefolds and 
fourfolds. In F-theory model building the base manifold B is a hypersurface in a four- 
dimensional toric ambient space. The fourfolds are complete intersections in a six-dimensional 
toric space. The associated lattice polytopes live in four- and six-dimensional integer lattices 
and typically have a large number of points. It is in general not possible to do calculations 
without computer support. There exist several software packages which are useful for par- 
ticular aspects in toric geometry. In this section we will mostly focus on the program PALP 
[29]. Before that, let us mention some other useful programs: Schubert by Katz and Str0mme 
is a Maple package for calculations in intersection theory. TOPCOM [35J computes triangu- 
lations of point configurations. Singular [39] is a powerful computer algebra program which 
is optimized for calculations with polynomial rings, such as the intersection ring. A recent 
addition is cohomCalg [22] which can compute line bundle-valued cohomology classes over 
toric varieties. 

Let us now discuss some features and applications of PALP [29] , which stands for "Package 
for Analyzing Lattice Polytopes". It consists of several programs. 

• poly . x computes the data of a lattice polytope and its dual if the polytope is reflexive. 
The input can be either a weight matrix or the points of a polytope in the M-lattice 
or the N-lattice. Apart from the polytope data poly.x computes Hodge numbers of 
the associated Calabi-Yau hyper surf aces, information about fibrations and other data, 
poly.x has been extended with several features that include information about the 
facets of the polytope, data of Fano varieties and conifold Calabi-Yaus. In [51\ [52] this 
extension of PALP has been used to find new Calabi-Yau manifolds with small h ' 
which are obtained from known Calabi-Yau threefolds via conifold transitions. The full 



12 



set of options in PALP can be obtained with poly.x -h and poly.x -x for extended 
options. 

• The program nef .x can be used for complete intersection Calabi-Yaus. It takes the 
same input as poly.x and computes the polytope data, nef partitions and Hodge num- 
bers as well as information about fibrations. There are several extended options which 
include most notably the data of the Gorenstein cones (cf. [53] for the definition and 
construction in toric geometry) in the M/iV-lattice. 

• cws . x creates weight systems and combined weight systems of polytopes of dimension 
to be specified in the input. 

• class. x classifies reflexive polytopes by searching for sub-polytopes of a Newton poly- 
tope associated to a combined weight system. 

Apart from recent applications in F-theory model building, which we will discuss in the next 
section, PALP has been used in many other contexts. A data base of Calabi-Yau threefolds 
has been generated by listing all 473 800 776 reflexive polyhedra in four dimensions [54] , In 
view of the landscape problem in string theory the statistics of the polytope data is also of 
interest [55J . Some of the most recent extensions of PALP which we will mention below have 
already been used in [561 Q3J El] • 

3.5 Application to F-theory GUTs 

In this section we make the connection to F-theory model building and discuss how the 
calculations discussed in section 12.21 can be carried out explicitly. The approach discussed 
here is used in [131 [24"] . Our aim is a systematic construction of a large class of examples of 
global F-theory models. The first step is the construction of the base manifold B. In [13] we 
have obtained a set of geometries by systematically constructing weight matrices associated to 
point and curve blowups on Fano hypersurfaces in P4. In [21] we have considerably extended 
this class of models by defining hypersurfaces in a subset of the toric ambient spaces described 
by the 473 800 776 reflexive polyhedra in four dimensions [51]. Concretely, we have restricted 
ourselves to configurations where the N-lattice polytopes have at most nine points. As one 
can check for example at [57], there are 1088 such polytopes. We used PALP to recover the 
toric data of the ambient space and considered all non-negatively curved hypersurfaces in 
these ambient spaces. In order to be able to perform the calculations outlined in section 12.21 
we must compute the intersection ring and the Mori cone. We have achieved this by using an 
extended version of poly.x [17]. The following additional features have been implemented: 
processing of non- Calabi-Yau hypersurfaces by specifying the hypersurface degrees as input 
parameters, a calculation of the maximal triangulations of the iV-lattice polytope, calculation 
of the Mori cone and the Stanley-Reisner ideal, and calculation of the intersection ring with 
the help of Singular. Using this data we can identify del Pezzo divisors, check the existence 
of a decoupling limit and compute the topological properties of matter curves and Yukawa 
points. In [21] we have analyzed a total number of 569 674 base manifolds. The resulting 
geometries are available at [55]. 

The next step in the calculation is to construct the Calabi-Yau fourfold which is an 
elliptic fibration over the base B. The toric data of X4 is obtained by extending the weight 
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matrix of B. Schematically, this looks as follows: 

3 2 1 

* * wn 

* * 

* * w m \ 

Here the Wij denote the entries of the weight matrix associated to B. The *-entries m the 
extended weight matrix have to be chosen in such a way that the fiber coordinates x, y are 
sections of Kg 2 and Kg 3 , respectively. These entries of the fourfold weight matrix contain the 
information about the hypersurface degrees of the base. Not every extended weight system 
will lead to a Calabi-Yau fourfold of the form ([1]). The calculations can be done using nef .x. 
Several problems can appear: first, there may be no nef partition and therefore our methods 
do not work. A second conceptual problem is that the polytope corresponding to the extended 
weight system is not always reflexive. Many of the combinatorial tools used in PALP are only 
valid for reflexive polytopes. Even though one might have a perfectly fine Calabi-Yau fourfold 
we cannot apply our technology to them. The third issue is of a technical nature: due to 
the complexity of the fourfold polytopes one may reach the software bounds of PALP which 
results in numerical overflows. For the 569 674 extended weight systems discussed in [23j we 
find only 27 345 reflexive fourfold polytopes which have at least one nef partition. Furthermore 
there are 18 632 reflexive polytopes without a nef partition, 381232 non-reflexive polytopes 
and 142 470 cases with numerical overflow. 

Having found a reflexive fourfold polytope with at least one nef partition is not enough 
to have a good global F-theory model. If we further demand that the base B has at least one 
del Pezzo divisor with a mathematical or physical decoupling limit the number of fourfolds 
decreases significantly. In addition we should also impose some constraints on the regularity 
of the base. Demanding that B is Cartier leaves us with 16 011 good models. Imposing the 
stronger criterion of base point freedom we are down to 7386 models. Focusing on these 7386 
good geometries we apply the constraint that the nef partition should be compatible with the 
elliptic fibration. This information can be extracted from the output of nef .x. This further 
reduces the number of geometries to 3978. 

Having found a good Calabi-Yau fourfold, we can construct a GUT model on every (del 
Pezzo) divisor. A toric description on how to impose a specific GUT group on a Tate model 
has been given in |6j. The Tate form ([2]) implies that the a n appear in the monomials which 
contain z n . We can isolate these monomials by identifying the vertex v z in (Vi,V2) that 
corresponds to the z-coordinate in the Tate model. All the monomials that contain z r are 
then in the following set: 

A r = {w k G A m : {v z ,w k ) + 1 = r} u z G V m , (26) 

where A m is the dual of V m , which denotes the polytope containing the z-vertex. The 
polynomials a r are then given by the following expressions: 

2 

a r - ^ c fc 11 11 Vi \ x =y=z=i (27) 

Now we can remove all the monomials in a r which do not satisfy the factorization constraints 
(UJ) of the Tate algorithm. In order to perform this calculation we have to identify the fiber 



(25) 
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coordinates (cc, y, z) and the GUT coordinate w within the weight matrix of the fourfold. We 
have applied this procedure to every del Pezzo divisor in the 3978 "good" fourfold geometries. 
Note that the procedure described above can destroy the reflexivity of the polytope, which 
happens in about 30% of the examples. For S'C/(5)-models we found 11275 distinct models^] 
with reflexive polyhedra, for SO(10) GUTs there are 10 832. L r (l)-restricted GUT models 
|14j can be engineered along the same lines. It turns out that [/(l)-restriction does not put 
any further constraints on the reflexivity of the polytope. 

4 Outlook 

In this article we have discussed how toric geometry can be used to construct a large number of 
geometries that can support global F-theory GUTs. Using this technology we could show that 
elementary consistency constraints greatly reduce the number of possible models. However, 
due to computational constraints, we did not quite succeed in systematically listing all possible 
F-theory models within a class of geometries. Such an endeavor would require substantial 
changes in the computer programs we are using. It is actually quite remarkable that we could 
make use of PALP for Calabi-Yau fourfolds and non-Calabi-Yau threefolds, since this goes 
beyond what it was originally designed for. 

Let us present a list of suggestions to extend PALP in order to improve the applicability 
to the current problems in mathematics and physics and to make it more accessible for users. 
The original purpose of PALP was to solve a classification problem for polytopes. Over the 
years it has been adjusted and extended in order to be applied to specific problems. Many of 
the basic routines that were implemented to tackle some special questions could be used in 
much more general contexts but cannot be easily accessed. Therefore a better modularization 
of the software is necessary in order to have flexible access to these basic routines. Another 
problem of PALP is that one has to specify several parameters and bounds such as the number 
of points in a polytope in a given dimension at the compilation of the program. It would be 
practical to have fully dynamical dimensions in order to work with a precision tailored to the 
problem at hand without recompiling. 

A fundamental change would be to step away from the description of polytopes and 
instead use the ray representation which has the full data of the cones. This is necessary if 
one wants to deal with non-reflexive polytopes. A further extension which has already been 
partially implemented is to include triangulations, intersection rings and even the calculation 
of Picard-Fuchs operators needed for mirror symmetry calculations into PALP. The ultimate 
goal is to have an efficient and versatile program which can be used for toric calculations of 
all kinds without having to rely on commercial software. Finally a detailed documentation of 
all the features of PALP would be helpful |59| . 

As for the search for F-theory models, an extended version of PALP would hopefully help 
to overcome the problems of non-reflexivity and overflows we have encountered in |24] . Apart 
from finding new examples for physics applications one might also attempt a partial classi- 
fication of Calabi-Yau fourfolds. Enumerating all toric Calabi-Yau fourfolds may be out of 
reach or even impossible but for finding all models of type ([I]) one can at least give a prescrip- 
tion for the construction: take each of the 473 800 776 reflexive polyhedra in four dimensions 
and put in all non-negatively curved hypersurfaces that are not Calabi-Yau. Then construct 

2 Since the procedure has been applied to all del Pezzo divisors in a given base geometry not all these models 
may have a decoupling limit. 
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fourfolds which are elliptic fibrations over these base manifolds. A rough estimate shows that 
this procedure would yield 0(lO n ) fourfold geometries. 
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